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Abstract 

A measurement of tlie ultrahigli-energy (UHE) cosmic neutrinos at a km^-size neutrino 
telescope will open a new window to constrain the 3x3 neutrino mixing matrix V and probe 
possible new physics. We point out that it is in principle possible to examine the non-unitarity 
of V, which is naturally expected in a class of seesaw models with one or more TeV-scale 
Majorana neutrinos, by using neutrino telescopes. Considering the UHE neutrinos produced 
from the decays of charged pions arising from pp and (or) collisions at a distant astrophysical 
source, we show that their flavor ratios at a terrestrial neutrino telescope may deviate from the 
democratic flavor distribution 0j : c;/)^ : 0j = 1 : 1 : 1 due to the seesaw-induced unitarity 
violation of V. Its effect can be as large as several percent and can serve for an illustration of 
how sensitive a neutrino telescope should be to this kind of new physics. 
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1 Introduction 



The solar [1], atmospheric [2], reactor [3] and accelerator [4] neutrino experiments have provided us 
with very convincing evidence that neutrinos are massive and lepton flavors are mixed. In the basis 
where the flavor eigenstates of charged leptons coincide with their mass eigenstates, the phenomenon 
of neutrino mixing can simply be described by a 3 x 3 unitary matrix V which links the neutrino 
flavor eigenstates {i/^, u^, v^) to the neutrino mass eigenstates {v^, v^^ v^: 

v,2 v^A\uA . (1) 

A full parametrization of V requires 3 rotation angles {6121 ^13, ^23) 3 phase angles (5, p, a) [5]: 

^ '^13*^12 '^13-^12 -^13^ *''\ 

^ ~ ~"^12'^23 ~ '^12=^23=^13^' +'^12*^23 ~ '5l2'^23'^13^' '^13'^23 "^M ' (2) 
V +512523 ~ <^12<^235l36* ~^12^2'i ~ 5]^2C235l3^* ^^^3023 / 

where s^^ = sin6^j, c^j = cosO^j (for ij = 12, 13, 23), and Pjy^ = Diag{l, e*^, e'"^} is the Majorana phase 
matrix irrelevant to neutrino oscillations. A global analysis of current experimental data [6] points 
to ~ and 6'23 ~ 7r/4, a noteworthy result which has motivated a number of authors to consider 
the yU-r permutation symmetry and its breaking mechanism for model building [7]. 

Now that neutrinos can oscillate from one flavor to another, it will be extremely interesting 
to detect the oscillatory phenomena of ultrahigh-energy (UHE) cosmic neutrinos produced from 
distant astrophysical sources. IceCube [8] , a km^-volume under-ice neutrino telescope, is now under 
construction at the South Pole and aims to observe the UHE neutrino oscillations. Together with 
the under-water neutrino telescopes in the Mediterranean Sea (ANTARES [9], NESTOR [10] and 
NEMO [11]), IceCube has the potential to shed light on the acceleration mechanism of UHE cosmic 
rays and to probe the intrinsic properties of cosmic neutrinos. An immediate consequence of neutrino 
oscillations is that the flavor composition of cosmic neutrinos to be observed at the telescopes must 
be difi^erent from that at the sources [12]. By measuring the cosmic neutrino flavor distribution, one 
can determine or constrain the mixing angles (^125^135^23) ^^e Dirac CP-violating phase (S). A 
lot of attention has recently been paid to this intriguing possibility [13] — [16]. 

We aim to investigate the oscillation of cosmic neutrinos produced from the decays of charged 
pions arising from energetic pp and (or) p'j collisions at a distant astrophysical source (e.g., active 
galactic nuclei or AGN). For such a most probable UHE neutrino source, its flavor composition is 

0, : 0^ : 0, = 1 : 2 : , (3) 

where 0^ = + (for a — e, /i, r) denotes the a-neutrino flux at the source. As the distances 
between the astrophysical sources and the terrestrial detectors are much longer than the typical length 
of solar or atmospheric neutrino oscillations, one may average the UHE cosmic neutrino oscillation 
probabilities and arrive at 

PaP^Pi^a^^p)=Y.\Vj'\V^^\'. (4) 

i=l 
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This result is also valid for the anti-neutrino oscillations; namely, Pap = Pii^a ~^ ^fs) — ^ap 
a, P — e, II and r. Therefore, the neutrino fluxes at the detector can be calculated from 

Given Eq. (3) together with the condition jV^J = \V^-\ (for i = 1,2,3) [17], it is easy to show that 
the flavor distribution of UHE cosmic neutrinos has a democratic pattern at neutrino telescopes: 

0j : 0j : 0^ = 1 : 1 : 1 . (6) 

Note that \V^j] — imphes either ^^^3 = and ^23 — '^/^ (CP invariance) or S — ±7r/2 and 
— n/ A (CP violation) in the standard parametrization of V as shown in Eq. (2). These two sets 
of interesting conditions can be realized from the so-called tri-bimaximal [18] and tetra-maximal [19] 
neutrino mixing scenarios, respectively. 

One has to bear in mind that (pj : (pj : cp'^ = 1 : 1 : I depends on two idealized hypotheses: 
the astrophysical source of UHE neutrinos satisfies 0g : 0^ : 0^ = 1 : 2 : and the 3x3 neutrino 
mixing matrix V satisfies \V^^\ = \V^^\. Previous works have extensively analyzed possible deviations 
from the democratic fiavor distribution of UHE cosmic neutrinos at neutrino telescopes by taking 
account of the energy dependence, uncertainties in the neutrino mixing angles, contaminations to 
the canonical production of i/g's (i^gS) and u^'s (z^^'s) from tt'^'s, and different sources of UHE cosmic 
neutrinos [13] — [16]. 

We shall concentrate on the standard pion-decay source of UHE neutrinos, whose flavor composi- 
tion has been given in Eq. (3), to explore the effects of non-unitarity of V on the flavor distribution 
of such cosmic neutrinos at a terrestrial neutrino telescope. This investigation is new and makes 
sense, because V is naturally expected to be non-unitary in a class of seesaw models with one or 
more TcV-scale right-handed Majorana neutrinos. We find that the democratic fiavor distribution in 
Eq. (6) can be broken at the percent level as a consequence of the unitarity violation of V. Although 
such a small effect is hard to be observed in any realistic experiments in the foreseeable future, it 
does illustrate how sensitive a neutrino telescope should be to this kind of new physics. 

2 Unitarity Violation at Neutrino Telescopes 

If the tiny masses of three known neutrinos {u^, U2, u^) arc attributed to the popular seesaw mechanism 
(either typc-1 [20] or type-11 [21]), in which there exist a few heavy (right-handed) Majorana neutrinos 
N^, then the 3x3 neutrino mixing matrix V must be no n- unitary. The effect of unitarity violation 
of V depends on the mass scale of N-, and it can be of O(10~^) if A''^ are at the TeV scale [22] — an 
energy frontier to be explored by the LHC. Indeed, a global analysis of current neutrino oscillation 
data and precision electroweak data yields some stringent constraints on the non-unitarity of V, 
but its effect is allowed to be of O(10~^) [23] and may have some novel implications on neutrino 
osciUations [24]— [26]. 

In the presence of small unitarity violation, we write the neutrino mixing matrix as V — AVq, 
where is a unitary matrix containing 3 rotation angles {9i2, ^13, ^23) 3 phase angles like that 



3 



given in Eq. (2), and A is a quasi- identity matrix which can in general be parametrized in terms 
of 9 rotation angles 9^^ and 9 phase angles (for i = 1,2,3 and j = 4,5,6) [24]. For simplicity, 
here we adopt the expression of A shown in Eq. (11) of Ref. [24] and take Vq to be the well-known 
tri-bimaximal mixing pattern [18] without any CP-violating phases. Then we obtain the non-unitary 
neutrino mixing matrix V = AV(^ as follows: 



/ 



2 



-j.{l-W, + 2X) 
\j={l-W,-2Y + Z) 



where 



for i — 1,2, 3; and 



1 



(1 - W,) 

X) 

W^ + Y + Z) 



'i4 



'i5 



'i6 



1 
V2 



{l-W^-Z). 



(7) 



(8) 



Y — S-y^S^^ + S15S35 + 516^36 , 
Z — 524-^34 + ■^25'^35 + '^26'^36 ■ 



(9) 



Here s. 



sinO^j and s^^ 



e^^i^s^j have been defined, and higher-order terms of s^j have been 
neglected. The mixing angles in can at most be of 0(0. 1), but the CP- violating phases 5^^ are 
entirely unrestricted. If both and are switched off, the tri-bimaximal neutrino mixing pattern 
will be reproduced from Eq. (7). With the help of Eqs. (4) and (7), we arrive at 



P.. = 



9-y^- 

l-\{W, + W,) + lReX, 

l-\{W, + W,)-'^-{ReY -2ReZ) , 

7 14 2 

w _ _ReX , 

18 9 ^ 9 

1 _ ! (ly^ + _ 1 (ReX - Rer + 2ReZ) , 
7 14 2 

--Wj + - (Rey - 2ReZ) . 



, , (10) 

18 9 ^ 9 ' ' ^ ' 

For the canonical astrophysical source of UHE neutrinos under consideration, we definitely have 
= {l/3,2/3,O}0Q, where 0q denotes the total initial flux. It is then easy to get the flavor 
distribution at a terrestrial neutrino telescope: 



K = — 



10 

3 
3 

00 

3 



2W^ 



4 

-ReX 
9 



^ {W^ + %W^) - ^ReX 

- {2W^ + + 9Wj) - ^ReX 
9 9 



(11) 
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The democratic flavor distribution of 0^ (for a — e, n, r) is clearly broken. Because of the non- 
unitarity of V, the total flux of UHE cosmic neutrinos at the telescope is not equal to that at the 
source: 



E 



;i2) 



This sum is apparently smaller than 0q, and it approximately amounts to O.960Q if Wi ~ 0.01 (for 
i — 1, 2, 3). Some comments are in order. 

(1) Note that ReX receives the most stringent constraint from current experimental data, l-'^l < 
7.0 X 10~^ [23]. Hence the dominant effects of unitarity violation on (f^ come from W^. The breaking 
oi (j)^ : (fy^ : (pj — 1 : 1 : 1 can be as large as several percent. Although the strength of unitarity 
violation is very small and certainly difficult to be observed in realistic experiments, it does illustrate 
how sensitive a neutrino telescope should be to this kind of new physics. 

(2) Note also that the oscillation probabilities of UHE cosmic neutrinos are actually given by 
= P^js/ (^^^)Qa(^^^)/3/3 (for a,j3 = e,fi,T) in the non-unitary case, where the production 

of i/q, and the detection of are both governed by the charged- current interactions [23]. Given 
the canonical source of UHE neutrinos, i/g's are generated from the decay of muons, and thus the 
charged-current interaction involves two lepton flavors (i.e., e and /x). But i/^'s can be produced from 
two channels: one is the decay of charged pions and the other is the decay of muons. The former 
involves only one lepton flavor (i.e., /i). Hence one should take care of the normalization factors when 
doing specific calculations of the cosmic neutrino fluxes for a specific neutrino-telescope experiment. 
For the simple pattern of V taken above, the normalization factors can be explicitly written as 



/2Wi X* 

X Z* I . (13) 

\ Y Z 2W3. 




(3) The unitarity violation of V under discussion is ascribed to the existence of heavy Majorana 
neutrinos in seesaw models and usually referred to as the minimal unitarity violation [23]. In contrast, 
the existence of one or more light sterile neutrinos and their mixing with three active neutrinos may 
also violate the unitarity of V. Using S^j to denote the matrix elements of active-sterile neutrino 
mixing, we can express the averaged probabilities of UHE cosmic neutrino oscillations as 

Pap = P{^a - = E + E \Saj?\Sp,? , (M) 

i=l j=l 

where a and /3 run over e, /i and r, and 

EIK./ + El'5a/ = 1, (forn = l,2,---) (15) 

i=i j=i 

holds. Eq. (15) shows the apparent unitarity violation of V induced by hght sterile neutrinos. Two 
observations have been achieved in Ref. [27]: (a) for small active-sterile mixing (i.e., I^'^^l <^ 1), the 
effect of non-unitarity of V at neutrino telescopes is very small and quite similar to that obtained in 
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Eq. (10); (b) for large hitherto-unconstrained mixing between active and sterile neutrino species (i.e., 
^ 1)) the existence of light sterile neutrinos might significantly modify the democratic flavor 
distribution of UHE cosmic neutrinos at neutrino telescopes. At present, however, we have to admit 
that there is no strong experimental or theoretical motivation to introduce light sterile neutrinos into 
the standard model. 

For illustration, we simply assume that there is only one heavy Majorana neutrino, which can be 
accommodated in the minimal type-II seesaw model [28]. In this case, we are left with three mixing 
angles (^^^4, ^24; ^'34) and three CP-violating phases (^^^4, ^24, ^34) characterizing the unitarity violation 
of V. As done in Ref. [14], three working observables at a neutrino telescope can be defined: 



R 
R 



^1 



Re « 


I 1 r 


2~36 


R, « 


1 1 r 
i - H 

2 36 L 


Rr « 


1 1 r 
i - H 

2 36 L 



Rr = iT^Zj.T ■ ^^^"^ 

n + K 

In the unitarity limit where V takes the tri-bimaximal mixing pattern, one can easily obtain R^ = 
R^ = R^ = 1/2, a result which is equivalent to the democratic flavor distribution. With the help 
of Eqs. (8), (9) and (11), we are able to evaluate the above flux ratios in the presence of unitarity 
violation: 

24s^4 — 15^24 — 9^34 — 12S]^4S24 COS Q 
12s^4 — 2IS24 + 9S34 — 6S;^4S24 cos Q 

12si4 + 6^24 - I8S34 - 6S14S24 cos g , (17) 

where g = 5^4 — 524 ^hc higher-order terms of s^j (for ij = 14, 24, 34) have been neglected. Taking 
into account the experimental constraints [23], we have numerically calculated the allowed regions of 
these working observables in Figure 1, where the phase angle g varies freely in the range g E [0, 27i]. 
Two comments are in order: 

• The deviation of R^ (for a = e,fi,T) from its value in the unitarity limit (i.e., R^ = 1/2) is at 
most at the 0.1% level. There are two obvious reasons for this result: (a) there exist significant 
cancellations among the contributions of three mixing angles to the flavor ratios; (b) the mixing 
angles s^^^ and S24 are strictly constrained by |X| = S;^4S24 < 7.0 x 10^^. 

• In more general cases with two or three heavy Majorana neutrinos, the above constraint can 
be loosened. Taking two TeV-scale Majorana neutrinos for example, we can obtain s^j ~ 0.1 
(for i = 1, 2, 3 and j = 4, 5) when the destructive interference between S14S24 •Si5'S25 terms 
takes place in X (see Eq. (9) and switch off the contribution of S16S26 ■^)- 

While a neutrino telescope is expected to identify different flavors of UHE cosmic neutrinos, it is also 
expected to measure the total flux as precisely as possible. A notable feature of unitarity violation 



6 



of V is that the total flux at the detector is not equal to that at the source, and such a discrepancy 
may be as large as several percent shown in Eq. (12). 



3 Comments on cosmic neutrino decays 

So far we have assumed cosmic neutrinos to be stable particles and studied their flavor distribution at 
neutrino telescopes. Now let us make some comments on cosmic neutrino decays and their possible 
signatures at neutrino telescopes. It is actually not unnatural to speculate that massive neutrinos 
are unstable and can decay into lighter neutrinos and other massless particles. If neutrino masses 
arise from spontaneous breaking of the global {B — L) symmetry, for example, then ~^ + X 
decays may take place, where x is a Goldstone particle (i.e., Majoron) [29]. A more exotic scenario, 
in which massive neutrinos may decay into unparticles, has also been proposed [30] . 

Here we consider a rather simple case: the decay products of UHE cosmic neutrinos are invisible, 
implying that the initial neutrinos simply disappear. When the neutrino source spectrum falls with 
energy in a sufficiently deep way, the daughter neutrino will also have negligible contributions to the 
total neutrino flux. Then the resultant neutrino flavor distribution at neutrino telescopes is simply 
given by [16, 31] 

0j:0j:0^ = |KiM^^ir:|V;ir, (18) 

provided is the lightest neutrino mass eigenstate (and thus stable). Note that Eq. (18) holds in 
the assumption that the heavier neutrinos V2 and completely decay into and invisible (massless) 
particles. If the neutrino mixing matrix V is not unitary, as illustrated in Eq. (7), then the flavor 
distribution at neutrino telescopes reads 

(/)J : 0j : 0^ = 4 (1 - 2W^) : (1 - + 4ReX) : (1 - 21^3 - 4Rey + 2ReZ) . (19) 

It is straightforward to compute the flavor ratios deflned in Eq. (16). In the unitarity limit, we have 
— 2 and i?^ — R^ — 1/5; and in the non-unitary case with only one heavy Majorana neutrino, 
we obtain 



R 



e 



" 5 + 25 

1 1 

R^ 



2Sl4 - S24 ~ "534 + 4S14S24 COS Q — 4S14S34 COS ■& + 2S24S34 COs(^ - "(9) 
1 1 



4si4 - 5^24 + S34 + 2OS14S24 COS g + 4S14S34 COS t!) - 2S24S34 cos(^ - •??) , 

r- . „^ 4s^4 + S24 - 5^34 - 4Si4S24COS^> - 20Si4S34COS1? + 10S24S34COS(£> - 'I?) , (20) 

5 25 L J 

where ^ = ^^^4 — ^24, ^9 = Si^^ — ^34, and higher-order terms of s^^- have been neglected. The allowed 
regions of three flavor ratios are plotted in Figure 2, where the phase angles g and vary freely in 
the range [0, 27r] . Two comments are in order: 

• Different from the case discussed in section 2, here the deviation of R^ from its value in the 
unitarity limit (i.e., R^ — 2) can be as large as 4%. In comparison, the deviation of R^ or R^ 
from its value in the unitarity hmit (i.e., R^ — R^ — 0.2) can be at the 0.1% level. 
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• It is worth mentioning that additional terms involving KeY and KcZ are present in Eq. (19), 
compared to Eq. (11). On the other hand, since s^^ or is confined to a very small value, the 
non-unitary CP-violating phase g can hardly affect the flavor ratios in Eq. (17). In the decay 
scenario, however, both the phases g and -d can significantly contribute to R^. 

We see that the flavor distribution of UHE cosmic neutrinos in the decay scenario is quite different 

from that in the standard neutrino oscillation picture. In particular, the democratic flavor distribu- 
tion of UHE cosmic neutrinos at neutrino telescopes is badly broken even if the condition | V^J = \ V^^\ 
(for i — 1,2, 3) is satisfied. 

4 Summary 

Assuming that UHE cosmic neutrinos are produced from the decays of charged pions arising from 
energetic pp and (or) p'j collisions at a distant astrophysical source, one may expect a democratic 
flavor distribution 0j : 0^ : 0!;"^ = 1 : 1 : 1 at neutrino telescopes if either ^^^^3 = and ^^23 = 7r/4 (CP 
invariance) or 6 = ±7r/2 and = 7r/4 (CP violation) are satisfied in the standard paranietrization 
of V. A lot of attention has been focused on small perturbations to the above conditions such that 
the resultant flavor distribution is no more democratic. We have explored a novel possibility, in 
which V is non-unitary and its non-unitarity is induced by heavy Majorana neutrinos as expected 
in a class of TeV-scale seesaw models, to examine the flavor distribution of UHE cosmic neutrinos 
at a terrestrial neutrino telescope. We have shown that the effect of unitarity violation on the flavor 
ratios 0j : 0^ : (f)J can be as large as several percent. We have also made some brief comments on 
cosmic neutrino decays and illustrated the relevant flavor distributions at neutrino telescopes. 

A measurement of the flavor distribution of UHE cosmic neutrinos is certainly a big challenge to 
IceCube and other neutrino telescopes. In the long run, however, we hope that neutrino telescopes 
can play an interesting role complementary to the terrestrial neutrino oscillation experiments in 
understanding the intrinsic properties of massive neutrinos and probing possible new physics. 

This work was supported in part by the National Natural Science Foundation of China. 
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Figure 1: Allowed regions of the flavor ratios (7?^, and {R^, R^), where the density of points is 
generated by scanning the possible ranges of s^j (for ij = 14, 24, 34) according to a flat random number 
distribution (i.e., s^j G [0,0.1] and S14S24 < '''■0 ^ 10~^ based on current experimental constraints on 
the non-unitarity of V). 
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Figure 2: Allowed regions of the flavor ratios (R^, R^) and {R^, R^) in the neutrino decay scenario, 
where the density of points is generated by scanning the possible ranges of s^j (for ij = 14, 24, 34) 
according to a fiat random number distribution (i.e., s.^^ G [0,0.1] and S14S24 < 7.0 x 10~^ based on 
current experimental constraints on the non-unitarity of V) . 
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